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Abstract 

Ground state of the quantum Jahn- Teller model with broken ro- 
tational symmetry was investigated by the variational approach in 
two cases: a lattice and a local ones. Both cases differ by the way 
of accounting for the nonlinearity hidden in the reflection-symmetric 
Hamiltonian. In spite of that the ground state energy in both cases 
shows the same features: There appear two regions of model parame- 
ters governing the ground state: the region of dominating selftrapping 
modified by the quantum effects and the region of dominating phonon- 
assisted tunneling (antiselftrapping) . In the local case (i) the effect of 
quantum fluctuations and anharmonicity due to the two-mode correla- 
tions is up to two orders larger than contributions due to the reflection 
effects of two-center wave function; (ii) the variational results for the 
ground state energy were compared with exact numerical results. The 
coincidence is the better the more far away from the transition region 
at the E®e symmetry where the variational approach fails. 
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Introduction 



Jahn- Teller model as a protopype model for phonon removal of degeneracy 
of electron levels in complex molecules [1] was investigated mainly in its 
rotational symmetric E®e local version with electron coupling to two degen- 
erate intramolecular phonon modes, one antisymmetric and one symmetric 
against reflect ion. Importance of focusing on the lattice version of the model 
increased due to JT-based structural phase transitions in some recently dis- 
covered high-T c superconductors and manganese-based perovskites [2], [3]. 

For the reflection symmetric two-level electron-phonon models with linear 
coupling to one phonon mode (exciton, dimer) Shore et al. [4] introduced 
variational wave function in a form of linear combination of the harmonic 
oscillator wave functions related with two levels. Two asymmetric minima of 
effective polaron potential turn coupled by a variational parameter respecting 
its anharmonism by assuming two-center variational phonon wave function. 
This approach was shown to yield the lowest ground state energy for the 
two- level models [4], [5]. 

The peculiarities due to reflection phenomena are also expected in the 
case of linear coupling with two phonon modes. In contrast to the model 
with coupling to one phonon mode, the oscillations mediated by the anti- 
symmetric mode are nonlinearly coupled to the symmetric phonon mode. 
Therefore, in order to improve the harmonic oscillator variational ansatz for 
this model, it was necessary to introduce additional variational parameters 
(VP) considering correlation of both phonon modes. This was performed by 
Lo [6] for the E®e JT model and by Lo et al. [7]) for the dimer model. 
The correlation was found most effective in the region of parameters where 
the competing localization (polaron) energy and the derealization transfer 
energy (tunneling) were comparable. 

In crystals exhibiting high space anisotropy the rotational symmetry of 
Jahn- Teller molecules can be broken. Therefore it is reasonable to investigate 
JT model generalized by breaking the rotational symmetry (while saving the 
reflection symmetry): we assume different coupling strengths for the onsite 
(intralevel) and interlevel electron-phonon couplings. Such a model can also 
be considered as a generalization of the exciton-phonon or the dimer-phonon 
model by assuming the electron tunneling phonon assisted. JT model with 
the broken rotational symmetry will involve the coupling of two minima as 
well as the coupling of the symmetric and antisymmetric phonon modes. In 



the lattice model, however, in contrast to the local case, the mode correlation 
appears to be of a marginal importance (Section III). As a consequence, our 
ansatz for the variational wave function of the local model will involve (i) 
the reflection VP introduced by the construction of the two-center phonon 
wave function, and (ii) correlation VP respecting coupling of the symmetric 
and antisymmetric phonon modes. The question arises of the importance of 
these VPs and respective effects regarding their relevance in different regions 
of the model parameters. Formulation of the variational ansatz of the local 
case is presented in the Section IV. Calculation of the ground state energy of 
the lattice model and related discussions are the contents of the Sect. III. and 
V. Reliability of different variational alternatives was checked by comparison 
with results of exact numerical simulations for the local model. 

2 Extended (lattice) generalized Jahn- Teller 
model 

We investigate ID lattice of spinless double degenerated electron states lin- 
early coupled to two intramolecular phonon modes described by Hamiltonian 

H = n (&lAn + V2) + £(c^l, n + M^n 

n,i=l,2 n 

-/?(&£,„ + b 2 , n )a xn ) - | E (Rij + R -^n ■ (1) 

Z n,j=l,2 

where i — 1,2 are phonon annihilation operators, and the Pauli matrices 
o\ n represent two- level electron system. They satisfy identities [ai n , <Jj n ] = 
i&kn, I = x,y,z, representing 1/2-pseudo-spins related to the electron densi- 
ties in a usual way, i.e. <J xn = l(Ci jn C 2 ,n + C2,n c l,n), °yn = j- {A,n c 2,n ~ 4,n c l,n) , 

e Z n = i(<i, n ci, n - 4,n c 2,n), h = \{A^\,n + 4,n c 2,n) is a unit matrix, and 
Cj „ are electron annihilation operators. The operator R±ij = e ±ipa of the 
displacement by a lattice constant ±a acts in both the electron and phonon 

Space, R±l,jfn = fn±lR±l,j- 

In terms of the creation-annihilation electron and phonon operators the 
Hamiltonian can be cast as follows: 

H = YJ 1 2 ^in h in + h + ?(ct,„Ci,„ - c\ n C 2 ,n){b\ n + b ln ) 
n i=l,2 Z A 

-^(c\ n c 2n + c\ n c ln ){b\ n + b 2n ) - tJ E ( c L c i,«+i + H - c -)] ■ ( 2 ) 

A A 3=1,2 



For P = —a, the interaction part of (1) 



/ b\ n + bi n , bin + ^2n 
Kbln + hn, -(b\n + hn) 



) 



(3) 



yields the rotationally symmetric E®e form [1] with a pair (an antisymmetric 
and a symmetric under reflection) of double degenerated vibrations. This is, 
e.g., the case of Cu ++ ions with d 9 configurations in high-T c cuprates [1],[2]. 

Taking a ^ (3 removes the degeneration of the vibronic states breaking 
the rotational symmetry of the electron-phonon interactions, the model still 
staying within the class of JT models [1], [2]. 

The dispersionless optical phonon mode b\ splits the degenerated unper- 
turbed electron level (j = 1, 2) while the mode 6 2 mediates the electron 
transitions between the levels. This latter term represents phonon-assisted 
tunneling, a mechanism of the nonclassical (nonadiabatic) nature as well as 
is the pure tunneling in related exciton and dimer models. 



Evidently, Hamiltonian (1) (a ^ (3) is reflection-symmetric, G^G'f^H = 



phonon 1 is antisymmetric under the reflection, phonon 2 remains symmetric. 

In addition, the transfer part of (2) exhibits SU (2) symmetry of the left- 
and right-moving electrons (holes). 

Let us note that the quantum phonon assistance of the electron tunneling 
(/3-term in (2) and (1)) constitutes the difference of the model from the related 
dimer and exciton quantum models where instead of (3^2{b\ n + b2 n )(Jxn of (1) 



there stands Aa xn , where A is the distance between the levels [4], [5]. 

The local part of (1) can be diagonalized in electron subspace by the 
Fulton-Gouterman unitary operator [8] U n = U2, n Ui, n -> where 



where Gffl 




n 




(5) 



as follows 




+ ^ EK & ln + Kn) ~ P{b\ n + b2n)G 1 , n )]I n . 



(6) 



n 



On the other hand, in the transfer term 



T 



(7) 



n 



there appears a nondiagonality 



Ki,±l — [(1 + Gi n Gi n ±i)I n + (1 — Gi n Gi ! n±i)G2 n O'zn\ 
X [(1 + G 2n G 2n ±i)I n + (1 — G^nG^nilVrrn] • 



(8) 




-i 



a x , and 



The diagonal terms of (8) represent the polaron transfer within one level 
while the off-diagonal ones represent the interlevel polaron transfer through 
the lattice. Evidently, the contribution of the off-diagonal terms proportional 
to 1 — Gi n Gi. ri +i is much smaller when compared with those proportional to 



Because of nonconservation of the number of coherent phonons, they are 
able even in the ground state to assist electron transitions between the levels. 
In the Hamiltonian (6), the operator G\ n = (— l) 6 i« 6ln (5), highly nonlinear 
in the phonon-1 appears mediated by phonons 2. It introduces multiple 
electron oscillations between the split levels mediated by continuous virtual 
absorption and emission of the phonons 1. The effect is analogous to Rabi 
oscillations in quantum optics due to photons [9]. Let us note that Rabi 
oscillations assist both the interlevel onsite and intersite electron transitions 
mediated by the electron transfer T. 

Further we shall investigate by variational approach the two variants of 
the model (1) - the lattice and the local one (considering T = 0). 

Shore et al [4] and Wagner et al [5] proved for exciton or dimer models 
coupled to one phonon mode that the two-center wave function generalized 
to an asymmetric nonunitary ansatz with a variational parameter r] of a form 



(where C is a normalization constant), yields lower ground state [4], [5]. 
So, rather generally we define phonon wave functions 4>^ 



1 + Gi n G 



¥ p) = -^={l + r]pG)\l)<p {p \ 



(9) 



= D l (±C l )D 2 (C 2 )S 2 (r 2 )S 1 (r l )S l2 (±\)\0). 



(10) 



Here, the generators of variational displacements are defined 



D i ( V )=exp[Q(bl-b l )], 



(11) 



and those of squeezing 

S t (r t )=exp[r t (bl 2 -b*)}, (12) 

for i — 1, 2. In Eq. (6) coupling of the phonon modes 1 and 2 occurs, 
therefore one includes into the ansatz (10) also the correlation generator 

S 12 (\) = exp[\(b\bl-b 1 b 2 )]. (13) 

with a correlation VP A. 

The form (10) is written here seemingly for the local case, while the 
most general form of it should contain the dependencies of all variational 
parameters on n (site number in the lattice): 71 — > 7i m (n) etc, showing the 
displacement of the mode 1 at the site m due to the electron at the site n, 
Ijmiji) = ^Eljqin) exp(-iqma). We take 

ljq(n) = 7y= exp{iqna) -> 7 jm (n) = 7j<W, (14) 

where 7, are independent of n. Eq. (14) indicates that the phonon dis- 
placement accompanies the electron at the site n. The same is valid for the 
squeezing, r m (n) = rb m ^ n and mixing parameter A. When we consider the 
ground state of the lattice model, we omit the electron and phonon dynamic 
terms by setting k = 0, q = 0, thus even in the lattice case we are left with 
the effective local molecule Hamiltonian, the only difference being the addi- 
tional contribution of the transfer terms (those containing T). As it is shown 
later, in the lattice model we can set A = without worsening considerably 
the variational results, while for the local case T = introducing this pa- 
rameter considerably improves the results. In what follows the ground state 
energy as function of optimized variational parameters r/, 71,72,?", A will be 
determined. 



3 Ground state of the lattice model 

First we consider the lattice model setting A = 0. 

By averaging H = H L + H T ((6), (7)) over the phonon wave functions (9) 
with (10) one obtains for the site Hamiltonian (6) the local part in the form 

= H ph + H int , (15) 



N 
where 

fi / , , m o 1 + V 2 - 2ve~ 8ri exp(-2|7!| 2 ) . |2r 

H ph = - cosh4r + 1 / + n ' ' FV ' 1 ' 7l 2 7 + 

2 1 + i] 2 + 2i] exp(— 2|7i| 2 ) 

+^|7 2 | 2 /, (16) 



TT _ rr , „ ^(l-?7 2 ) „ 
-Hint = tt a + Hf3= — ^—r / ,~ | 2 \ 7l^ - 

/5[(l+7 ? 2 )exp(-2|7i g | 2 )+2^ /j 



1 + r] 2 + 2r7exp(-2|7i| 2 ) 

where H a and Hp are cc- and /3-dependent parts of the interaction Hi nt (17). 
From the transfer Hamiltonian Ht (7) there remains 

^ = -Texp(-^)M = 2(# T1 7 + # r2 ^ + # T3 <7 2 + H TA a x ). (18) 
where M: 



M = 



1 + Ef)I + (1 - Ei)E 2 a zn . (1 + £ 2 )/ + (1 - Ei)<r x , (19) 



and Ei are given by 

Ex = exp (-2|7!| 2 ) , £ 2 ee exp (-2| 72 | 2 ) . (20) 

Expressions (20) in the ground state are independent of q and n because 
of the form of 7i g (n) = 7, exp(— igna) (14). This substantially simplifies 
subsequent calculations leaving and Ci as functions of 7? independent on 
n. The transfer terms Jf^, i = 1,2,3,4 in (18) with (19) are expressed as 

Hti = ~{1 + El){\ + £ 2 ), tf r2 = -|(1 - £ 2 )(1 - E 2 2 )E 2 , 

H T3 = ~(l-Et)(l + Ei)E 2 , H TA = -^(1 + Ef)(l-El). (21) 

Here, Hti,Ht3 and Ht 2 ,Ht4 are diagonal and off-diagonal terms of (18), 
respectively. 

The effective polaron potential in (15) (with (16)- (20)) as a highly non- 
linear function of ^1 and 7 2 is visualized in Fig. 1. 

The potential exhibits two sets of minima related to two competing 
ground states: 

(i) Two nonequivalent broad minima related to both the levels (10) at ±71 7^ 
and 7 2 close to 0; 

(ii) One narrow minimum at 71 close to 0, 7 2 > 0, where both levels approach 
close together. This minimum develops at growing (3; evidently its behaviour 
depends also on T value because of nonlinearity of the Debye- Waller factor 
M ((19) and (20)). 



Figure 1: The effective potential for x — P/ot = 0.5, /i = 2 and T — 1,Q — 
0.05 as function of 71, 72. 



The ground state energies related to these two sets of competing minima 
were calculated numerically. The result of the numerical minimalization of 
the diagonalized form of energy (15) 

r 1 V 2 

E = H ph + Hp + H Tl - (H a + H T3 ) +H T4 -H T2 (22) 

can be written symbolically as 

Eg(p, X, M/T) = E({-y lq , 7 2g , r, r]} min )/T. (23) 

Here, the index min denotes the optimized values. The model parameters 

« 2 (3 Q 

^=2^' X= a' T (24) 

are parameters of the effective interaction, asymmetry and nonadiabaticity, 
respectively. Energy Eq (23) is in scaled units, renormalized by T. The re- 
sults of the numerical evaluations of the ground state energy (23) are depicted 
in Figs. 2a, 2b different model parameters. 

One can distinguish there two regions depending on the value of x with 
different behaviour of ground state in each of them: 

(i) The ground state pertaining to the lower broad minimum at 71 < and 
72 pa with a small reflection part at 71 > is referred to as a "heavy" 
region. It corresponds to predominantly intralevel "heavy" polaron which 
is represented by a two-peak wave function, both peaks representing a har- 
monic oscillator (~ exp[— (x=F7i) 2 ]) displaced by the value ±71 (as it is seen 



from the form of the ansatz (10)). The "heavy" polaron is dominant at % < 1 
where the broad minimum 71 < is dominating. 

(ii) Close to x = 1, the energies of two minima go very close together (their 
difference is of the order of the phonon energy Q). They drop to one narrow 
minimum which represents a new ground state. Close to x — 1? continuous 
transition to a new ground state occurs. It stabilizes at x > 1> where (op- 
timized) 71 value is close to and 72 > 0. This region is referred to as a 
"light" polaron region because, owing to the abrupt decrease of 71 ps at 
X > 1, the effective mass of the intralevel polaron drops to almost its free- 
electron value, i.e., the polaron selflocalization vanishes in the "light" region. 
Because of this "undressing" , the transport characteristics of the excited elec- 
tron would increase. Moreover, due to the tiny distance between the levels, 
their coupling takes place by the exchange of virtual phonons 1. At suitable 
conditions, in the excited state, when both levels are occupied by electrons of 
opposite spins, the mechanism of virtual phonon exchange implies the pairing 
of electrons, i.e. formation of "light" bipolarons . 




Figure 2: The ground state energy (28) in the x~^ plane at ji = 2.5 (a) and 
in the plane x ~ A* — 1) (b)- The seltrapping dominated phase at x < 1 
and the tunneling dominated phase at x > 1. 

The ground state energy, especially its behavior dependent on pairs of 
parameters x> & and X> A 4 i s illustrated in details in Fig. 2. While being 
weakly fi-dependent, the energy strongly decreases with x inside the "light" 
phase (Fig. 2a). The position of the phase line is slightly shifted from x — 1 
at n — to higher values of x with increasing Q. This is consistent with the 
fact that the phonon fluctuations are most effective when the difference of 
the energies of the phases is of the order of the phonon energy. 

The ground state energy in the "heavy" region (x < 1) is independent of 
X, its decrease inside the "light" region (x > 1) is dependent on the effective 
coupling n (Fig. 2b). The energy decrease due to fi is stronger in the "light" 
phase (depending on x) than in the "heavy" phase. 



Figure 3: The range of relevance of the reflection parameter 77: Difference 
AEg = Ea(r) = 0, r = 0) — Eg, for % = 1, Eq is the exact ground state (23). 
For x < I5 the difference increases with x reaching its maximum at x — 1- 
At x > 1 h drops to zero, i.e. the narrow minimum "light" phase is resistant 
against 77. 

The region of importance of reflection measure 77 and r is illustrated in the 
Fig. 3. There we show the difference between the ground state energy with r 
and i] omitted and ground state with four variational parameters calculated 
numerically. However, contribution due to r E G (i],r = 0) — E G ~ 2.10 -3 
of the maximum value in Fig. 3. The said parameters turn out to be most 
relevant in E®e JT case (x= 1)- 

The effects of r\ and r on the ground state is apparent only at moderate 
couplings in the "heavy" region, reaching their maximum at x — 1- The 
narrow minimum (the "light" region) is resistant against rj. 

In order to demonstrate important properties of displacements (/i and 
Q dependence) it is sufficient, except of the region of importance of the 
reflection (Fig. 3) inside the "heavy" region apparent for 72, to calculate 
them in the limit 77 = 0. We shall use Eq. (22) approximated for the cases 
of both "heavy" and "light" polaron and obtain implicit expressions for 7^ 
which however are good to visualize their behaviour in both regions: 

(i) "heavy" polaron (72 small, E 2 ~ l,H T2 ,H Ti ~ 0): 



72 is small except of the region of fluctuations visualized in Fig. 4. One can 
see, that the selflocalization due to phonons 2 in the "heavy" region ~ —7! 
clearly implies the "cave" in the ground state energy due to the reflection 
measure i] (anharmonicity) of the ground state (Fig. 3). 
(ii) "light" polaron (7! small, E 1 l,H T2 ,H Ta ps 0): 

— ae 4r _ 3 

* * 2 ^ (l + ^ + ^) - 10 '- * ~ 2773(2 ^ (26) 

In the " light" region the fluctuations of 71 are missing as well as the fluctua- 
tions of the energy. This is consistent with the above result of the resistance 
of the narrow minimum against rj. 

For both "heavy" and "light" polaron a dependence of 7$ on the nonadi- 
abaticity parameter Q/T appears. It implies the dependence of the Debye- 
Waller factor and consequently of the polaron mass on the phonon frequency 
Q. This can be thought of as an analogy of the isotope effect at zero temper- 
ature. 

In the "heavy" phase, electron transitions mediated by phonons 2 to the 
upper level enhance fluctuations ~ 72 which mix phonons 2 with phonons 1 
and contribute to the fluctuations of the ground state energy of the heavy 
region. This is the reason for the similarity of the results in Fig. 3. 



4 Ground state of the local model 

Now we consider the same Hamiltonian, but taken for the local molecule (T = 
0). From the outset we include the additional variational mixing parameter 
A into the Ansatz. The understanding of the importance of introducing the 
mixing parameter in the local model in the contrast to the lattice one can 
be gained if one examines the expressions for Hamiltonians, in particular, 
those terms which contain the mixing of two phonon modes. The mixing is 
contained mostly in the "local" term with /3, 

/3Q 2 expi7r(bfb 1 ). 

The transfer term T also contains the mixing, but it always enters the 
expression in the form T(1+E 1 E 2 ) : thus for the transfer term its contribution 
is always shadowed by the larger term ~ 1. Thus, if T 7^ and large enough, 
its contribution is always dominant over mixing term f3. But if T = 0, we 
are left with the only nonlinear coupling term 3Q2S and accurate accounting 
for mode mixing becomes crucial. 

The variational mean value of the part of the Hamiltonian (6) in the state 
(9) renormalized by Q yields (Appendix A) 



(H) 1 

l—L = _(cosh4ri + cosh4r 2 ) cosh2A + (Ci + Cf) 
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exp 



2C( 



cosh 2A 



{— tanh 2A sinh 2A cosh 2(r x + r 2 ) cosh 2(r x — r 2 ) 



C cosh 2A 

+Ci[(exp2(ri +r 2 ) - exp(-2(ri + r 2 )) cosh4A)(l + tanh 2 2A) 
+2exp(— 2(ri + r 2 )) sinh 4A tanh 2A] cosh2(ri — r 2 ) 

+ %— (sinh4n - sinh4r 2 ) - 2Ci(tanh(2A)e 2 ( r2 - ri )C 2 + Ci) 
cosh 2A 



(27) 



2a 



2£ 



(1 +r/ 2 )[C 2 - Citanh2Aexp2(r 2 - n)]- 



exp 



2C? 



cosh 2A 



cosh 2A 



+ 277C2 



where 



ex ( 2C * ) 

Ci = Ciexp(-2ri), C = 1 + 2ry eXF _, ro s , + r/ 2 . 



cosh2A 1 ■' • (28) 
From Eq. (28) one can see that A causes correlations of the selftrapping 
and tunneling dominated regions. 
For A = 0, Eq. (28) turns to 



^ A ~° = |(cosh4n + cosh4r 2 ) + ^ (l - 2r/e" 8ri exp(-2C 2 ) + rf) £ 

2/3 



2 

«Cl + ^(i-^ 2 )Ci - 



nc 



c 

[(l + 7 7 2 )exp(-2C 1 2 ) + 27 7 ]C2 



where 



C = l + 2r/exp(-2C 1 2 )+r/ 2 . 
For r] = 0, C = 1 and (28) becomes 

(-H")t?=0 



(29) 
(30) 



n 



i(cosh 4ri + cosh 4r 2 ) cosh 2A + Ci + C2 + 77C1 

2Cl 



cosh 2A 



2/3 exp 
--(C2-Citanh2Aexp2(r 2 -r 1 ))— cogh2A 



(31) 



The ground states of (28), (29) and (31) were found by minimalization of the 
expressions against the involved VPs. The respective ground state energies 
E opt , £^= , E°^1q and -E^=o,A=o,o P t will be compared mutually and with the 



exact value from numerical simulation in order to find out importance of the 
variational parameters rj, A in different regions of model parameters x = 
(reflection) and fi = (effective e-ph coupling). The place of E®e JT model 
will clearly come out as an important special case. 

5 Discussion of the numerical results 

In lattice electron-phonon models the nonadiabaticity parameter is the ra- 
tio of phonon frequency and the band width, p. The energy is scaled by T 
and the effective coupling parameter (ratio of the polaron energy and of the 
phonon frequency) is \xt = a 2 /2QT). In the present local model, the nonadi- 
abaticity parameter is the ratio of the frequency and the coupling parameter 
Vt/a. The energy was scaled by Q and the ratio of polaron energy and the 
frequency is then a 2 /2Q 2 . Therefore, the reduction of the ground state en- 
ergy of the local model due to selflocalization is stronger than that of the 
lattice model by the factor [il/ = T/Q > 1. The energy of the tunneling 
(3 2 /Q 2 , is comparable with the polaron energy at (3 2 /fl 2 = /jll. Because close 
to x — 1 there is £i ~ 0, the effects of all variational parameters for x > 1 
except of (2 vanish. Quantum fluctuations due to r 1; r 2 and anharmonicity 
due to A are concentrated in the crossing region close to /x = 1 and x = 1- 




Figure 4: The numerical ground state wave functions at fi — 2 and x = 0-9 
(a) X = 1 (b) and x = 1-5 (c). 



The ground state in the phase plane x an d shown in Fig.4 exhibits 
two phases separated by the crossover line close to % = 1. It means, that 
the effective polaron potential exhibits two minima which reflect positions 
of two levels governed by the model parameters /i and x- The minima co- 
incide within the border of the regions close to x — 1- The phase x < 1 
is x independent, selftrapping dominated, with quantum fluctuations due to 
A,77,ri,r 2 . The phase x > 1 is the phonon-2-assisted tunneling dominated 
region with continuum emission and absorption of virtual phonons-1. The 
phonon exchange couples the levels within one minimum. The minimum is 
much more sensitive to the change of model parameters /i, x as well as to 
quantum fluctuations r±,r2 and A. 

Electron in the selftrapping dominated region is trapped by the phonons-1 
but due to the interactions mediated by phonons-2 the electron can fluctuate 
to the higher level. Due to the reflection symmetry of the phonons-2 contin- 
uum oscillations of the electrons at simultaneous emission and absorption of 
phonons-1 occurs. These oscillations couple the levels and so the electrons 
into pairs. This mechanism was described in a recent paper [10] for a lattice 
model. 

From (31), for small Xi one can approximately take (2 ~ 0, Ci ~ —fi and 
for A we are left with the result A ~ /ixexp(— 2/x). This dependence can be 
recognized in Fig. 5. 




Figure 5: The difference of E(r] = 0, \ opt ) - E(r] opU \ opt ). 



Comparing Figs. 6a and 6b one can see that the mode correlation induced 
anharmonicity A (Fig. 6a) is by one order larger than the contribution to 




Figure 6: The difference of E(i] = 0, X op t) — E(i] opt , X op t) including the mode 
correlation (a), without the correlation (b). Strong suppression of the reflec- 
tion effects by the mode correlation at /j < 1 is evident. The correlations 
flatten the barrier between the zones and the role of the parameter r\ is 
minorized (Note the difference of the scales of Figs, (a) and (b)). 



the anharmonicity of the reflection level mixing 77 Fig. 6b. The correlation 
is most effective for weak effective couplings fi at x ~ 1- F° r large /i it 
contributes only close to x = 1, where it reveals maximum for all /x In 
order to check the validity of the calculations using variational approach we 
performed also numerical diagonalization of the Hamiltonian in the phonon- 
1, 2 space. We limit ourselves in N\ phonon-1 states and N2 phonon-2 states, 
thus state vector is Ni x N 2 dimensional. As numerical diagonalization results 
show, about 20-50 phonon states are sufficient for convergence. We show the 
results of numerical diagonalization of Hamiltonian matrix as function of x 
for /i — 1 and /i — 4. In the first case we took 20x20 state vector, while in the 
latter case to achieve satisfactory convergence (especially for the tunneling- 
dominated region when % > 1) we had to increase the number of phonons-2 
up to 50. 



6 Quantum fluctuations in the local model 



Fluctuations of phonon-1 coordinate AQ\ = (Qf) — (Q1) 2 , Qi = ^(b\ + h) 

and the conjugate momentum P 1 = =j=(b\ — hi) for 77 = can be easily 
calculated analytically. The wave function of the phonon-1 has the Gaussian 
form 

4>[x) ~ exp —7——- — — r 

YK ' F [ cosh (2 A) exp(4ri) 

In fact, it is squeezed displaced harmonic oscillator if one does not take 
into account the parameter r\\ Introducing A does not invoke higher order 
nonlinearities with respect to phonon-1 coordinates. Thus (APi) = 0, like it 



should be for an harmonic oscillator, and the expressions for second momenta 
have the form: 

AQl = - exp 4ri cosh 2A, AP X 2 = - exp(-4ri) cosh 2A (32) 

and 

(AP 2 AQ 2 ) 1/2 = I cos h 2A. (33) 

The shape of fluctuations due to the correlation A (33) closely follows the 
shape of the curve in Fig. 7, having its maximum close to x — 1 an d /i ~ 1. 



\J < A P 2 A Q 2 > 




Figure 7: The product of uncertainties of the coordinate and the conjugated 
momentum of the phonon-1 for rj = 0. 

The similarity of the product of uncertainties of phonon coordinate and 
momentum displayed in Fig. 7. illustrates anharmonic source of the contribu- 
tion due to the mode correlation. The correlations A manifest themselves in 
(32) and (33) as a phonon anharmonicity. It is noteworthy to emphasize its 
quantum origin (compare the note below Eq. (2)). The problem with rj ^ 
needs special consideration and will be considered elsewhere. 

7 Conclusion 

A suitable choice of the variational wave functions for various electron- 
phonon two-level systems is a long-standing problem in solid state physics 



as well as in quantum optics. For two-level reflection symmetric systems 
with intralevel electron-phonon interaction the approach with a variational 
two-center squeezed coherent phonon wave function was found to yield the 
lowest ground state energy. The two-center wave function was constructed 
as a linear combination of the phonon wave functions related to both levels 
introducing new VP. 

This symmetry implies coupling of the levels mediated by continuous vir- 
tual emission and absorption of phonons accompanying the electron tunneling 
between the levels. This is an analogy of Rabi oscillations in quantum optics. 
Investigation of referring properties (tunneling rates, Debye- Waller factors) 
was performed mostly for the models with onsite (Holstein) electron-phonon 
interactions in dimer or exciton-phonon models and also for E®e JT model. 
The methods used there were based on either combination of variational 
approach and unitary transformations or numerical diagonalization. 

In two-level models with phonon assisted tunneling there appears coupling 
of both phonon modes mediated by Rabi oscillations. Therefore, also the 
interlevel quantum correlations of the phonons must be taken into account. 
This brings into the variational ansatz additional quantum VP. 

The two-level system is characterized by strongly nonlinear effective po- 
tential provided by phonons. The nonlinearity is governed by the bare Hamil- 
tonian parameters exhibiting a crossover from an asymmetric double poten- 
tial well with two broad nonequivalent minima pertaining to the levels to a 
regime of one narrow minimum when the broad minima coincide into one 
potential well. This effect is due to the effective attraction of the levels by 
virtual exchange of phonons. Much effort was expended also to improve 
one-center variational approach by including correlations of the concerning 
two squeezed coherent phonon modes of the JT model. When compared to 
the lattice case, the local version of the model provides qualitatively similar 
effective potential. This conclusion was to be expected: in the symmetric 
case minima of the nonlinear effective potential related with two levels are 
close together and so the quantum fluctuations there are most effective. We 
considered both the mixing and squeezing correlation VP in order to find re- 
gion of relevance of both parameters as functions of the interaction constants. 

The support from the Grant Agency of the Czech Republic of our project 
No. 202/01/1450 is highly acknowledged. We thank also the grant agency 
VEGA (No. 2/7174/20) for partial support. 

Appendix A 

We used following formulas [11] 

D 1 ( m )S 1 (r 1 ) = S 1 (r 1 )D 1 (f ll ), fj, = j^e"* (34) 



S l 1 (ri)biSi(ri) = 6icosh2ri + fr{sinh2r 1 (35) 
Sfa 1 (A) hSn (A) = h cosh A + b\ sinh A (36) 

(0|^(r)g}(^)ex P (A6t))= (coBh 1 2r)1/2 



x exp 



A 2 11. 
— tanh(2r) — Ar^(tanh2r — 1) — -rf + -(tanh2r — 



(37) 



Jin 

(0|5 1 (C)I>i(i7)6} ,B > = !^<0|Si(OA (^)exp(A6 t 1 ))| A=0 , (38) 
S 12 (A) = Tt S 1 (X/2)S 2 (-X/2)T (£) (39) 

r (i)a:)^(i)^C:;t)' 

T f |0) =T|0) = |0) (41) 
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